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0> ' Abstract 

Mh, Using a three-dimensional hydrodynamic simulation of the collision and an equa- 

qj ' tion of state containing a first order phase transition to the quark-gluon plasma, we 

study thermal photon production for Au + Au collisions at Eiab = 11-5 AGeV and 

^^ . for Pb + Pb collisions at 160 AGeV. We obtain surprisingly high rates of thermal 

H ; 

Cu ^ photons even at the lower energy, suggesting that, contrary to what was expected so 

far, photon production may be an interesting topic for experimental search also at 
the Alternating Gradient Synchrotron. When applied to the reaction S + Au at 200 
AGeV, our model can reproduce preliminary data obtained by the WA80 Collabora- 
tion without having to postulate the existence of an extremely long-lived mixed phase 
as was recently proposed. 
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In ultrarelativistic heavy-ion collisions one of the most challenging goals is to study 
a possible phase transition between nuclear matter and the quark-gluon-plasma (QGP), 
consisting of deconfined quarks and gluons. This expectation is supported by lattice 
gauge calculations and phenomenological models. Photons are a very promising probe 
in the experimental search for the QGP which is expected to exist for a brief period of 
a few fm/c before the majority of the final state hadrons is emitted. Particles which 
interact only electromagnetically may be regarded as free particles, carrying information 
on the early stage of the collision. In order to describe the space-time development of 
a nucleus-nucleus collision, very often hydrodynamical models are used, which are based 
on the concept of local thermal equilibrium. For heavy ion collisions at the Brookhaven 
Alternating Gradient Synchrotron (AGS) and at the CERN Super Proton Synchrotron 
(SPS) this condition may be fullfilled because of the high degree of stopping. Thus, high 
compression of nuclear matter has to be considered and energy densities near the critical 
value for the deconfinement phase transition may be reached. 

In the present paper, we study thermal photon production in a realistic hydrodynamical 
framework and make predictions for measurements at AGS and SPS energies. We focus 
on the reactions Au + Au at Eiab = 11.5 AGeV and Pb + Pb at 160 AGeV. We also apply 
our model to S -|- Au collisions at 200 AGeV and compare our results to preliminary data 
on photon production obtained by the WA80 Collaboration ||l|. 

At the SPS, lead beams will be available at the end of 1994 and direct photon pro- 
duction measurements will be performed by the WA98 Collaboration 0]. Furthermore, it 
is intended [^ to test the PHENIX-detector Q, which is designed to detect photons with 
high accuracy and resolution, at the AGS in 1996. This means that our predictions for 
Pb + Pb at 160 AGeV and for Au + Au at 11.5 AGeV can be checked in the near future. 

The photon production rates in thermalized matter have been calculated both for a 
QGP and for a hadron gas. In a QGP of temperature T, the rates for a baryochemical 



potential /x = can be expressed as [^: 

where i^-y and A: are the energy and three-momentum of the photon, and where a = 1/137 
and as are the electromagnetic and the temperature dependent strong coupling constant, 
respectively. From lattice results as{T) can be parametrized as [^: 

"'^^^ " (33-2n/)ln(8r/rc)' ^^^ 

where Tc is the critical temperature and nj is the number of flavors (below, we take 
Uf = 2). 

For a hadronic resonance gas at temperature T and baryochemical potential fj. = 0, the 
thermal photon rates were calculated in ||^. The authors of Q parametrized their results 
in the form 
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where Ej and T are measured in units of GeV. Eq. (|3|) contains the contributions from 
^i-resonances which were shown to be important in R and which had not been taken 
into account in earlier calculations of the rates |^ . 
In the mixed phase the rates read 
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where w{e) is the fraction of QGP at the energy density e. 

In order to obtain the single-inclusive spectra it is necessary to integrate the photon 
rates over the space-time region defined by the space-time evolution of the hot and dense 
matter, 

where the temperature field T{x) and the four velocity field u^{x) are obtained from 
a hydrodynamic description of the expanding matter. On the rhs of eqs. (|l|) and (0) 



the energy E^ has to be replaced by k^u^{x). In the framework of hydrodynamics the 
ultrarelativistic heavy ion colhsion can be divided into three stages: (i) the colhsion or 
compression stage, (ii) the expansion stage and (iii) the freeze-out stage. The space-time 
development of dense and hot matter in local thermal equilibrium is described by the 
equations of relativistic hydrodynamics: 

a^T^^ = 0, 5^(6n^) = 0, (6) 

where T^^ is the energy-momentum tensor, b the baryon- number density and u^ the four- 
velocity of the fluid element. Here we consider the case of an ideal fluid, i.e., we neglect 
dissipative effects. In this case, T^^^ = {e + P)u^u^ — g^'^P, where e is the energy density 
P the pressure. In order to solve the equations, one needs to specify the equation of state 
(EOS), which may be given, e.g., in the form P = P{e,b). The equations of motion are 
solved exactly and fully 3-dimensionally |^ with the numerical code HYLANDER. 

Our hydrodynamical simulation describes also the initial compression (or collision) 
stage where the propagation of shock waves which heat up the system is calculated fully 
3-dimensionally. The applicability of one-fluid hydrodynamics at this early stage is sup- 
ported by the fact that in heavy ion collisions at AGS and at SPS energies the mean free 
paths are small compared to the size and lifetime of the system and a considerable amount 
of stopping is expected^. 

We use an EOS given by a parametrization |Q of lattice-QCD results [^] where the 
hadronic phase is treated as a resonance gas. This EOS exhibits a first order phase transi- 
tion between hadronic matter and the QGP at a critical temperaturepl Tc = 200 MeV; it 
^RQMD simulations for central Au + Au collisions have found |10[| complete stopping up to (and even 

somewhat beyond) SPS energies. 

^The results in [|_1[ were obtained for a purely gluonic system. New results which take into account the 

effects of dynamical quarks yield lower values for the critical temperature, Tc ~ 150 MeV (see |12| and 

refs. therein). There is no EOS published yet for this case. 



corresponds to a baryocheinical potential /i = 0. Apart from the EOS our hydrodynamic 
model does not contain any free parameters. 

Fig. 1 shows the first 8 fm/c of the space-time evolution of central Au + Au collisions 
at Eiab = 11.5 AGeV (left column) and of central Ph + Pb collisions at 160 AGeV, in time 
steps At = 2 fm/c. Contour plots for the energy-density calculated with HYLANDER 
are displayed. Four characteristic regimes are considered. The energy density e = 0.14 
GeV/fm'^ corresponds to the interface between hot and cold nuclear matter, e = 0.25 
GeV/fm^ to the freeze out regime, the region 2.5 < e < 5.5 GeV/fm^ to the mixed phase 
and e > 5.5 GeV/fm'^ to the pure QGP phase, respectively. As can be seen in the figure, for 
both reactions a lump of pure QGP and a mixed phase with a lifetime of about 7 fm/c are 
produced. The total lifetimes of the thermalized matter are ~ 12 fm/c for Au + Au at 11.5 
AGeV and ~ 16 fm/c for Pb + Ph at 160 AGeV. Figs. 2a and 2b show the thermal photon 
spectra for these two reactions. The rates are surprisingly large even at AGS energies. To 
understand the origin of these results, we have plotted separately the contributions from 
the pure QGP phase, the mixed phase and the purely hadronic phase. For Au + Au at 11.6 
AGeV the contribution of the purely hadronic phase dominates the photon production, 
whereas for Pb + Pb at 160 AGeV it is comparable to the contribution of the pure QGP 
phase. The reason that the pure hadronic phase plays such an important role in photon 
production is twofold. Firstly, as can be seen from eqs. (|l]) and (|^, at temperatures 
T ~ 0.2 GeV the hadron gas outshines the QGP by a factor of about two. Secondly, the 
hadronic space-time volume exceeds the volumes occupied by the mixed phase and by the 
QGP. The fact that at SPS energies the QGP contribution can compete with that of pure 
hadronic phase is a consequence of the high initial temperatures (Tj ^ 300 MeV) at these 
energies. 

Preliminary data on photon production for S + Au collisions at SPS energies were 
■^At present there exist no reliable lattice data concerning the EOS at /x 7^ 0. 



recently presented by the WA80 Collaboration [|l|. The fact that the measured direct 
photon rates are surprisingly high has led to speculations concerning the existence of an 
extremely long-lived mixed phase with lifetimes of about 30 — 40 fm/c[0]. We have applied 
our hydrodynamic collision simulation to this reaction^ and plotted the resultant thermal 
photon rates in Fig. 2c. We find a remarkably good agreement with the WA80 data ||l|, 
without having to adjust any parameters or having to assume an extremely long-lived 



mixed phase as in Ref. [13|. In particular, it turns out that the hadronic rather than the 
mixed phase dominates photon emission. The space-time evolution for S + Au collisions 
is illustrated in Fig. 3. Due to the asymetry of the collision, the central rapidity region is 
shifted, which is of importance for the calculation of the photon spectra. 

So far we have not discussed the effects of baryon stopping on thermal photon pro- 
duction. At AGS and SPS energies, one expects nonvanishing values of the baryochemical 
potential due to the high baryon number densities in the central rapidity region. In Ref. 
[ p!7[| , photon emission rates from a QGP of temperature T and baryochemical potential 
H were calculated by means of the Braaten-Pisarski technique. The authors give the 
expression 

At fixed baryon number density and energy density, we have used a bag model EOS with 
Tc{fj, = 0) = 0.2 GeV to estimate the decrease in temperature and the resultant decrease 
in the photon rates [18, ^ near the critical curve of the deconfinement phase transition. 



Using eqs. (|^ and (^ we obtain reduction factors of about 2 both for the QGP and for 
*For the system S + Au one expects a smaller degree of stopping than for Pb + Pb or Au + Au. On the 
other hand, fits to rapidity and transverse momentum spectra of hadrons produced in 5 + 5 at 200 AGeV 
already indicate a considerable amount of stopping |14| which should be even more pronounced for S + Au. 
This suggests that our hydrodynamic collision simulation is also applicable for the latter reaction. Note, 
however, that recently the WA80 data have also been described M, uM under the assumption of Bjorken 
initial conditions. 



the hadronic component^. 

From the above considerations we conclude that the main result of this letter - namely, 
that there are experimentally observable rates of thermal photons to be expected both at 
the AGS and at the SPS - remains unaffected if one takes into account a finite baryochem- 
ical potential. 
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^ A reduction by a factor of ~ 2 for S + Au at 200 AGeV would imply that we underestimate the 
WA80 data. There are three possible effects which would lead to an increase of the photon rates and 
thus may explain the missing factor. An additional contribution is also expected from the hadronic 
processes Ai -^ -K-y, 61 -^ ■k-k'^^, and Ki -^ K-y [H. Furthermore at nonzero baryochemical potential the 
presence of baryons may open additional channels for thermal photon production. Thirdly, one expects a 
decrease of the rho meson mass as the temperature approaches the critical temperature of the chiral phase 
transition teoll. 



Figure Captions 

Fig. 1 Energy density contour plots in the (x, r) plane, for Au + Au collisions at Eiab = 
11.5 AGeV (left column) and Pb + Pb collisions at Eiab = 160 AGeV (right column). 
X is the coordinate in longitudinal direction, r the radial coordinate in the trans- 
verse plane. The plots were obtained in the equal velocity frame by applying the 
hydrodynamic collision simulation. 

Fig. 2 Single-inclusive photon spectra as a function of the transverse momentum kj_ for 
mid-rapidity photons, (a) for Au+Au collisions at Eiab = 11.5 AGeV, (b) for Pb+Pb 
at 160 AGeV and (c) for S + Au at 200 AGeV. The contributions from the pure 
QGP phase, the mixed phase and the pure hadronic phase are shown separately. 
The data points are from ref. ||l|. 

Fig. 3 Energy density contour plots as in Fig.l for S + Au collisions at Eiab = 200 AGeV. 
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This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411336vl 
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This figure "figl-3.png" is available in "png" format from: 
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